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Abstract 

In many circumstances the perfect fluid conservation equations 
can be directly integrated to give a geometric-thermodynamic equa- 
tion: typically that the lapse N is the reciprocal of the enthalpy h, 
(N = 1/h). This result is aesthetically appealing as it depends only 
on the fluid conservation equations and does not depend on specific 
field equations such as Einstein's. Here the form of the geometric- 
thermodynamic equation is derived subject to spherical symmetry and 
also for the shift-free ADM formalism. There at least three appli- 
cations of the geometric-thermodynamic equation, the most impor- 
tant being to the notion of asymptotic flatness and hence to space- 
time exterior to a star. For asymptotic flatness one wants h — > and 
N — ► I simultaneously, but this is incompatible with the geometric- 
thermodynamic equation. A first shot at modeling spacetime exterior 
to a star is to choose an idealized geometric configuration and then 
seek a vacuum-Einstein solution. Now the assumption of a vacuum is 
an approximation, in any physical case there will be both fields and 
fluids present. Here it is shown that the requirement that the star 
is isolated and the presence of a fluid are incompatible in most isen- 
tropic cases, the pressure free case being an exception. Thus there is 
the following dilemma: either an astrophysical system cannot be iso- 
lated or the exterior fluid must depend explicitly on the entropy. That 
a system cannot be isolated is another way of stating Mach's princi- 
ple. The absence of asymptotically flat solutions depends on: i)thc 
equation of state, ii)the admissibility of vector fields, and iii)the re- 
quirement that the perfect fluid permeates the whole spacetime. The 
result is robust against different choices of geometry and field equa- 
tions because it just depends on the fluid conservation equations and 
the ability to introduce a suitable preferred vector field. For example 
with spherical symmetry there is the preferred vector field tangent to 
the 3-sphere; furthermore for asymptotically flat spacetimes there is 
the preferred vector field tangent to the 3-sphere at infinity. The prop- 
erties of spherically symmetric geodesies are studied both to examine 
whether a fluid spheres geodesies can explain solar system dynamics 
and also to model hypothetical galactic halo with spherically symmet- 
ric fluids so as to produce constant galactic rotation curves. The rate 
of decay of fields and fluids is discussed in particular whether there 
are any non-asymptotically flat decays and what gravitation modifies 
the Yukawa potential to be. The Tolman-Ehrcnfest relation follows 
immediately from N — 1/h. The equations relating the enthalpy to 
the lapse have a consequence for the cosmic censorship hypothesis this 
is briefly discussed. 
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1 Introduction. 



1.1 A Relation Between Fluids and Geometry. 

The perfect fluid stress can be covariantly differentiated to give the perfect 
fluid conservation equations. In many cases these differential equations can 
be directly integrated to give a geometric-thermodynamic equation, which 
typically equates the Eisenhart (1924) Q Synge (1937) [§] fluid index w 
to the reciprocal of the lapse N. The Eisenhart-Synge index is essentially 
the fluids zero temperature enthalpy. In section ^ the index is calculated 
for several equations of state. The a— equation of state is a 2— parameter 
equation of state which describes polytropes. The /3— equation of state is 
a 1— parameter equation of state obtained from the a— equation of state 
by assuming the second law of thermodynamics for an adiabatic process, 
Tooper (1965) Q, Zeldovich and Novikov (1971) ||]. It gives the 7— equation 
of state in all cases except for 7 = where the pressure free (p = 0) case is 
not recovered; but rather /i = pln(-fe). For equations of state see also Eligier 
et al (1986) ||, and Ehlers ||. The index is not defined in the pressure free 
case, thus solutions such as the Tolman (1934) 0-Bondi (1947) || solution 
are not covered by description in terms of the fluid index. 

1.2 It's Application 

The main application of the geometric-thermodynamic equation is to the 
description of spacetime exterior to a star. It is shown in many cases that 
asymptotic flat solutions do not exist. This is taken to imply that the 
notion of asymptotic flatness as usually understood is physically simplistic. 
In the literature diagrams are constructed which are supposed to represent 
the causal spacetime of a collapsing star. These diagrams usually require 
that the spacetime is asymptotically flat, but the inclusion of a non- vacuum 
stress is often sufficient for this requirement no longer to hold. An example 
of how these diagrams can be qualitatively altered by infinitesimal matter is 
given by solutions to the scalar-Einstein equations which often have no event 
horizons, Roberts (1985) ||. To the lowest approximation the spacetime 
exterior to a star has no stress: the star exists in a vacuum. In order to 
take account of the matter that surrounds a star it is necessary to find an 
approximate stress which has contributions from planets, dust etc. . . There 
seems to be no systematic way of producing a stress which approximates 
such diverse forms of matter. When relativity is applied to macroscopic 
matter the stress is usually taken to be a perfect fluid, so that this is taken 
to be the form of the first order correction to the vacuum. Specifically the 
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stress is taken to be a spherically symmetric perfect fluid with 7— equation 
of state and the result generalized where possible. The nature of the surface 
of the star is left open as boundary conditions to interior solutions are not 
discussed. The assumed equations of state are essentially of one variable 
so that the pressure p does not depend on the entropy s, i.e. p = p(fi) 
not p = s) or in other words they are isentropic. Stars radiate and 
the radiation possess entropy whether there are equations of state that can 
describe this, and if so whether they are susceptible to a similar analysis 
as to that given here is left open: however it would be unusual for at the 
simplest level there to be asymptotically flat spacetimes, at the next level 
of complexity for there to be none, and at the full level of complexity for 
asymptotically flat spacetime to reappear. 



1.3 Asymptotic Flatness. 

It will be shown that many spacetimes with a perfect fluid stress do not 
have asymptotically flat solutions. Throughout it is assumed that the fluid 
permeates the whole spacetime and that the spacetime is of infinite extent. 
Also throughout it is assumed that a simple limiting process is appropriate: 
specifically as the luminosity radial coordinate tends to infinity the space- 
time becomes Minkowskian. This does not always happen, two examples 
are: Krasihski's (1983) Ji~0[ analysis of the Stephani universe where r = 00 
labels a second center of symmetry, and type B Bekenstein conformal scalar 
extensions of ordinary scalar field solutions, Agnese and LaCamera (1985) 
[ |11| |, and Roberts (1996) J12|| , which have bizarre asymptotic properties. The 
result follows from the conservation equations so that it is explicitly inde- 
pendent of the gravitational field equations used. The conservation equa- 
tions use a Christoffel symbol or a generalization of this. The connection 
depends on the metric which in turn can be thought of as a solution to grav- 
itational field equations. In this implicit sense the result can be thought of 
as depending on field equations. It might not hold if there are other fields 
coupled to the fluid. Similarly asymptotically flat solutions are rare for 
theories with quadratic Lagrangians, Buchdahl (1973) |l3| . The absence of 
asymptotically flat solutions might have application to the "missing mass" 



problem, see Roberts (1991) [14] and references therein. Some properties 
of exact perfect fluid solutions have been discussed by Delgaty and Lake 
(1998) @. In Bradley et al (1999) ||] it is shown that the Wahlquist per- 
fect fluid spacetime cannot be smoothly joined to an exterior asymptotically 
flat vacuum region. Boundary conditions for isolated horizons are described 



in Ashtekar et al (1999) |17]]. Models of stars in general relativity have been 
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discussed by Herrera et al (1984) @, Stergioulas (1998) @, and Nilsson 
and Uggla (2000) @. 



1.4 Sectional Contents. 

Section || introduces the stress, conservation equations, and the relationship 
between the enthalpy h and the Eisenhart-Synge fluid index cu. In section ^ 
it is shown that there are no asymptotically flat static fluid spheres unless 
the fluid index u — > 1 at infinity, and using Einstein's field equations there 
are no asymptotically flat static fluid spheres with 7— equation of state. In 
the non-static case for 7— equation of state there are no asymptotically flat 
solutions provided that 7/0,1 and certain conditions hold on the met- 
ric. For both static and non-static cases there might be asymptotically flat 
solutions for a— polytropes. In section |] it is shown for static spacetimes 
admitting non-rotating vector U a = (N, 0) and having 7— equation of state 
that the lapse N is inversely proportional to the fluid index N = 1/lo. For 

the non-static case, subject to N = and 7(ln(^ g^)),i = 0, the equation 

relating the lapse to the fluid index is u> = N^ 1 g^ 2 = ^~^T . These 
results can be used to show that there are no asymptotically flat fluid fill- 
ing spacetimes admitting the vector U a = (N, 0) with 7— equation of state, 
again also subject to certain conditions on the metric. The introduction 
of the vector U a = (N, 0) assumes that the fluid is non-rotating and that 
the spacetime admits a global time coordinate, unlike the vacuum Einstein 
equations, see for example Cantor et al (1976) 21], and Witt (1986) |p2| . In 



section [| a case against asymptotic flatness is presented. Outer solar system 
observations of orbital irregularities are discussed. Non-asymptoticness on 
length scales greater than the solar system, such as galaxies, is mentioned. 
The time-like geodesies for an arbitrary Newtonian potential are calculated. 
Modeling hypothetical galactic halos of "dark matter" with spherically sym- 
metric fluid solutions so as to produce constant galactic rotation curves is 
attempted. The the rate of decay of various fields are discussed. It is ar- 
gued that most perfect fluid spheres and some conformal scalar spheres rate 
of decay is in fact an increase prohibiting asymptotic flatness. There is 
the possibility of experimentally testing gravitational theory by measuring 
the deviation of the Yukawa potential from what would be expected in the 
absence of gravitation; how this might be done is briefly discussed, the pos- 
sibility of an actual test seems remote. Various onion models of spacetime 
surrounding the Sun are discussed. It is argued that non-asymptoticness 
implies that a system cannot be gravitationally isolated and that this sug- 
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gests a new formulation of Mach's principle: there are no flat regions 
of physical spacetime. The philosophy of what an "isolated system" 
entails is briefly discussed. In section || the Tolman-Ehrefest (1930) [p3|| 
relation is derived. Section [7| speculates on the relevance of the geometric- 
thermodynamic relation to cosmic censorship. 

2 The Enthalpy and the Eisenhart-Synge fluid In- 
dex. 

2.1 Perferct Fluids. 

The stress of a perfect fluid is given by 

T a/3 = (v + p) U a U/3 + p g af3 = nh U a Uj3+p g a/3 , U a U a = -1, (1) 

where /i is the fluid density, p is the pressure, n is the particle number, h 
is the enthalpy, and p + \x = nh. The unit timelike vector U a defines the 
geometric objects 

h a /3 = 9al3 + U a Uf3, U a = UatfljP , 8 = U" a , K af3 = U x -sh^h^, 

u a p = h*hpU[ x .s], cr aj3 = U( a .p) +U( a Uf3) - -9h aj3 , (2) 

called the projection tensor, the acceleration, the expansion, the second 
fundamental form, the rotation, and the shear, see for example page 83 of 



Hawking and Ellis |24|. The projection obeys U a h a p = and U a h a p = Up, 

also the acceleration obeys U a U a = 0. Formally the second fundamental 
form and its associated hypersurface only exist when the rotation vanishes. 
Transvecting the stress conservation equation rA with U a and h°L gives 
the first conservation equation 

" U a T?. p = HoJJ a + + p)U« a = /i + (fi + p)6 = (3) 

and the second conservation equation 

h a X.;P = (V + P)U« + h a % = 0, (4) 

respectively. These equations equate the derivatives of the vector field to 
the pressure and density. From a technical point of view, here we are inves- 
tigating when these equations can be integrated. It turns out that assuming 
a specific form of vector field - say hypersurface orthogonal U a = \(j)^ a is not 
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directly of much use, but rather assumptions about the form of the metric 
have to be made. The first law of thermodynamics can be taken in the 
infinitesimal form 

dp = n dh + nT ds, (5) 

where T is the temperature and s is the entropy. The Eisenhart[|l|-Synge|Q] 
fluid index is defined by 

after setting T = in || and integrating it is apparent that up to a constant 
factor at zero temperature oj = h. The index is also discussed on page 84 of 
Hawking and Ellis [pi|]. 

2.2 Polytropes 

The a— polytrope has equation of state 

p = a/ (7) 

and has 

dp = afifi^dfi, (8) 

or 

because of this the pressure is not an explicit function of two variables a 
and (5, but only one. The index and particle number corresponding to |7] are 

u = (l + a//" 1 )/^, n = + a^~ l )^P , (10) 

The (3— polytrope 0] has equation of state 

p = Kn~<, (11) 

where K is a constant and V = 1/n is the volume occupied by one bayron. 
For an adiabatic process (no exchange of heat) the second law of thermody- 
namics is 

dE , s 

>=-w <12) 

where E is the total energy density per unit mass E = fi/n. Then [T^ 
becomes 

P = n— , (13) 
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11 and ^ give 



P n~* = = (14) 

which in the case 7^1 can be integrated to give 

fi = —n\ (15) 
7-1 

where the constant of integration is taken to be zero. Using |TT], |l^ becomes 
the equation of state of 7— polytrope 

p=( 7 -l)At, 7/1, (16) 

which has index and particle number 

7-1 1 
to = /j 7 , n = 7/x^, 7^0,1, (17) 

In the pressure free case (7 = 1 in |l6[) the index ^ is not defined, an option 
is to replace p with (7 — l)fi in the definition || and then take 7 = 1 to 
obtain In (a;) = or to = 1, then the condition nu = fi + p gives n = [i. For 
the 7— equation of state the first || and second || conservation laws can be 
written in terms of fj,, where fi = w- 1 , and are 



/i + 7A i0 = O, (18) 

and 

lf iU a + (7 - l)h£np = 0, (19) 

respectively. The 7— equation of state has been derived under the assump- 
tion that 7/I. Perhaps the correct 7=1 equation of state for a (3— polytrope 
is found by putting 7 = 1 in 14 and integrating to give 

p_ 
K 



fi = pln[^); (20) 



however the speed of sound 



B " 'lnfel+1] , (21) 



dp V \Kj 

is 1 or the speed of light when p/K = 1, it is less than the speed of light 
for p/K > 1, and it diverges as p/K — > exp(— 1). That the speed of sound 
can take these values suggests that this equation of state is essentially non- 
relativistic. Some writers refer to ^ as dust, others call the pressure free 
case p = dust. ^ has index and particle number 

u> = (l + ln(|)J K , n = p(l + ln(£)J * . (22) 
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3 Asymptotically Flat Fluid Spheres 
3.1 Spherical Symmetry. 

The line element of a spherically symmetric spacetime can be put in the 
form 

ds 2 = -C dt 2 + Adr 2 + B d£ 2 . (23) 
Choosing the timelike vector field 

U a = (VC, 0,0,0), (24) 

the rotation vanishes and the projection tensor, acceleration, expansion, 
shear, and second fundamental form are 

K = h 6 = h* = i, 
u« = (o,^,o,o), 

1 , A B. 



VC"2A B' 

2A 2A 1 1 ■ A ■ 

arr = -^°» = -B*iM <r ** = Sjd(- A+ B B) > 

K rr = -i-U 

K ™ = lL K ^ = -\7c B ' (25) 

where the overdot denotes absolute derivative with respect to r as in U a = 
^r 2 -, but otherwise the overdot denote partial derivative with respect to 
time. Noting that dfi/dr = dt/dr fij = fi/y/C, the first conservation equa- 
tion H] becomes 

A-(M + P)(^ + |) = 0, (26) 
and the second conservation equation ||| becomes 

P' + (V + P)^ = 0, (27) 
only the r component is non-vanishing in the second equation. 
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3.2 The Static Case. 

In the static case the first conservation equation Eg vanishes identically and 



the second conservation equation 27 integrates to give 



" " 75' (28) 



the constant of integration is taken to be independent of and (j) and is 



absorbed into C, for example by redefining t. For the line element |23| to 
be asymptotically flat it is necessary that as r — ► oo, the line element ^ 
becomes Minkowski spacetime in other words as r increases C — * 1, A — * 1 



and B — > r 2 . Now from 28, C — » 1 implies that u; — * 1. Thus any static 
spherical fluid sphere with a well defined index not equal to or 1 cannot 
be asymptotically flat. To see this result in particular cases first consider 



the 7— equation of state. From 17 and 



3 



^ = 0^=^, (29) 

and as C — > 1, /x tends to a constant and thus the spacetime cannot be 
asymptotically flat; also the spacetime cannot be asymptotically DeSitter 
as this would necessitate \x tending to a constant time r 2 . In the pressure 
free case, the index is not defined and there are the asymptotically flat 
solutions given by Tolman M and Bondi M. Next consider the (3— equation 



of state, from 22 and 



C=(l + ln(|)J (30) 

now asymptotically as C — > 1, p — > K; however a constant value of p 
asymptotically is not consistent with asymptotic flatness, therefore there 
are no asymptotically flat solutions. Finally consider the a— equation of 



state, from 10 and 28 



C = (1 + a//" 1 )^, (31) 

in the case p — > 0, C — > 1 and there might be asymptotically flat a— polytropic 
spheres. The same results are obtained using the more general vector 

U a = (ay/C,yj(a?-1)A,0,0), (32) 

where a is a constant. 
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3.3 The Non-static Case. 



In the non-static case it is necessary to assume an equation of state in order 
to calculate a geometric-thermodynamic relation. The 7— equation of state 
is assumed. Then either from [l^ and |l^ and 25, or from 16 and 26 and ^7] 
the first and second conservation laws are 

M-7M^ + 5 J=0, (33) 

and 

( 7 _iy + 7/x g = , (34) 
respectively. The equation 

d[i = jldt + fi'dr 

- -#+f)-+^r£*- <35) 

can be integrated when 

(7 = 0, (36) 

and 

(AB 2 )' = 0, (37) 

to give 

H = uj^i = A+ifi+^C 5 ^), 7^0,1, (38) 

where the constant of integration have been taken to be independent of 9 
and cj> and is absorbed into the line element. The assumption (AB 2 )' = 
is coordinate dependent and holds rarely as for example (AB 2 )' ~ 4r 3 
for Minkowski spacetime in spherical coordinates whereas (AB 2 )' = for 
Minkowski spacetime in rectilinear spacetime. Taking the limits A, C — ► 1, 
B — > r 2 , for 7 > 0, /i — ► a constant, and for 7 < 0, fi diverges; thus there 
are no asymptotically flat solutions. The a— equation of state cannot 
be investigated without further information. Discussion of non-existence of 
time dependent fluid spheres can also be found in Mansouri (1977) [p6| . 

4 The Geometric-Thermodynamic equation in the 
ADM formalism. 
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4.1 Vanishing Shift ADM Formalism. 



In the ADM (-l,+3) [^7]] formalism with vanishing shift the metric is given 
by 



(39) 



where is the determinant of the 3— dimensional metric. The reason the 
shift is taken to vanish will become apparent later. The timelike unit vector 
field used here 



U a 

Ui; t 

U t -t 



-N.4 



u- - -— « (3) 



U, 



i:i 



(40) 



there are other choices such as U a = (— iV, 0), and also U a = (aN,bNi) for 
which the unit size condition U a U a = — 1 implies g^NiNj = a ^ l . For ^ 
the rotation vanishes and the remaining geometric objects |2| are 



hij — 9ij ■• 



N- 



N 



In {g 



(3) 



Ki. 



L «J ~~ x \?« ~~ 9ij,t- 

The first conservation equation |3] becomes 



am - (a* + p) ( ln \fg~^^J 



o, 



and the second conservation equation becomes 



Pi + ^ + P)^- = 0, 



(41) 



(42) 



(43) 



the t component of the second conservation equation 43 vanishes identically. 
If the shift is included in the above vector one finds 

2N 2 U® = + (N k N k ) ti + Ni(2N j;i - N k g {ikj} ), (44) 

and further calculation proves intractable. 
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4.2 Static Case. 



In the static case the first conservation equation vanishes identically and the 
second conservation equation integrates immediately and independently of 
the equation of state to give 

where the constant of integration has been absorbed into N. 



4.3 Non-static Case. 



In the non-static case assume the 7— equation of state has to be assumed to 
accommodate the first conservation law ||. With 7— equation of state 16 the 
conservation equations 42 and 43 integrate to give 



lo(3)|(7-l) 

N 



7/0,1, 



where in place of 36 and 37 



JV = 0, 



and 




0, 



(46) 



(47) 



(48) 



respectively. Constants of integration have been absorbed into the line ele- 
ment. Substituting the spherically symmetric values of the previous section 
times a function of sin 9 which has been taken to be absorb 
depend on the choice of velocity vector 



into 46 



gives 



able there. The equations |45| and |_ 
40, for example if a geodesic velocity vector is chosen then the acceleration 



vanishes and 45 and Hq do not hold. The conditions 47 and 48 do not appear 



to have an invariant formulation. There are three things to note. The first is 
that these derivatives do not occur in the covariant derivatives of the vector 
field EH and hence do not occur in the geometric objects 41. The second is 
that and 48 are satisfied if 



{**} = o, 



and 



(49) 
(50) 

The third is 

that |47| and |48| might solely be a gauge condition; but |47| puts on one 



{$*} = 0, 

respectively, as they only occur in these Christoffel symbols. 



constraint and 48 puts on three constraints totaling four, the usual number 
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of differential gauge constraints. The Plebanski-Ryten (1961) |25| gauge 
condition is 

(51) 



[(-9) w Ab = o, 



for w = \ this is the harmonic gauge condition. For a = t, 51 



is 



N 
AT3 



For a = x l , 51 



is 



0. 



0. 



(52) 



(53) 



For w ^ 0, |47| and [48| cannot be recovered, except for Minkowski spacetime 
in rectilinear coordinates. Thus the conditions ^ and 48 on the metric 
appear not to be an example of Plebanski-Ryten gauge conditions. It can 
be asked, is there a non-static geometric-thermodynamic relation which in- 
volves familiar gauge conditions instead of metric constraints such as [47] and 



48 



Inspection of [T^ and ^ with arbitrary vector field instead of |40| does 
not immediately give a choice of vector field for which application of the 
Plebanski-Ryten gauge [|l] simplifies matters enough for the problem to be 
tractable. 



4.4 7-equation of state and the ADM. 

For the 7— equation of state ^ becomes 

fi = N^, 7^0,1, (54) 

for the spacetime to be asymptotically flat the density \x must vanish asymp- 
totically implying that the lapse N must vanish, contradicting the assump- 
tion that the spacetime is asymptotically flat. For the 7— equation of state 



16, 45 becomes 

fi = N^g {3 ^ /2 , 7/0,1, (55) 

asymptotically /x — > r 2 and the spacetime cannot be asymptotically flat. For 
a— polytropes the static case ^ gives 

N=(l + ayP~ l ) 1-/3 , (56) 

and in this case it is possible for iV — ► 1 and \i — > simultaneously as r — > 00 
Thus for spacetimes where the rotation free vector 4C can be introduced, and 
subject to the caveats mentioned above for the non-static case: i)there are 
no asymptotically flat 7— polytropes except possibly for 7 = or 1, ii)there 
are no asymptotically flat fluid spacetimes unless the fluid index tends to a 
finite non- vanishing constant. 
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5 Against Asymptotic Flatness. 



5.1 Lenght Scales. 

On length scales from the outer solar system to cosmology there are observa- 
tions indicating that asymptotic flatness of the systems under consideration 
are not correct. It is known that the dynamics of the outer solar system 
have unexplained irregularities. For example from the figures of Seidelmann 
et al (1980) |E8| it appears that the irregularity in Pluto's orbit is that the 
RA increases by about 2 arcsec more than expected in 50 years, similarly 
the declination decreases by about 1 arcsec in 50 years. The irregularities 
are not neatly expressible by a single quantity, as for example the orbit of 
Mercury was prior to general relativity; but roughly this means that the or- 
bit is boosted by about 2 arcsec in 50 years. This makes the construction of 
theories to explain the irregularities difficult. In Roberts (1987) p9| the ef- 
fect of a non-zero cosmological constant was investigated in order to explain 
the irregularities of Pluto's orbit and it was found that the cosmological 
constant would have to be about 12 orders of magnitude bigger than the 
upper bound Zel'dovich (1968) |3t| finds from cosmological considerations. 
Axenides et al (2000) |3(| also discuss dynamical effects of the cosmological 
constant. 



Scheffer (2001) [31] and Anderson et al (2001) discuss dynamical irreg- 
ularities in the paths of spacecraft. The orbit of comets, Marsden (1985) 
p4| , Rickman (1988) pq ], and Sitarski (1994) |}(| have unexplained irregu- 
larities, for example at least 6 comets have unexplained forces acting toward 
the elliptic. Qualitatively this is exactly what would be expected from Kerr 



geodesies [37| page 363, 

In summary then, the bound and the marginally bound orbits 
must necessarily cross the equatorial plane and oscillate about 
it. 

but qualitatively the effect is many orders of magnitude out: on solar sys- 
tem length scales the Kerr modification of Schwarzschild geometry is intrin- 
sically short ranged. These solar system orbital problems might originate 
from the oblateness of the sun, Landgraf (1992) (j38|. There are theories 
which have gravitational potential with an exponential term and mass scale 
m p (mH/'mp) n , where mu is a typical hadron mass, mp is the Planck mass, 
and n = 0, 1, and sometimes 2. Satellite and geophysical data for n = 2 
theories show that they are not viable unless mu > W^GeV, Gibbons and 
Whiting (1981) [ |39[ . Other searches for an adjusted potential have been 
undertaken by Jarvis (1990) 
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5.2 The Exterior Schwarzschild Solution as a Model. 



The exterior Schwarzschild solution is a reasonable model of the solar sys- 
tem outside the sun. A fluid solution can be argued to be a better ap- 
proximation to the matter distribution as it takes some account of inter- 
planetary space not being a vacuum. Any exterior fluid spacetime would 
have different geodesies than the vacuum Schwarzschild solution, conse- 
quently the orbits of the planets would be different from that suggested by 
the Schwarzschild solution: how to calculate these geodesies for spherically 
symmetric spacetimes is shown below. The magnitude of the upper limit 
of the effective cosmological constant is about p\ = 10~ 16 g. cm. -3 , that 
this is too small to explain Pluto's irregular orbit was shown in Roberts 
(1987) [p9f| . Thus to explain Pluto's irregular orbit using a fluid the critical 
density must be larger than p\. p\ is much larger than the mean den- 
sity of interplanetary space which is of the order of 10 _29 g. cm. -3 (or 10~ 5 
protons cm." 3 ). The density of interplanetary matter is insignificant com- 
pared to the density contribution from the planets, for example for Jupiter 
AJupiter = i^Jupite r ry u 3 pitcr » 2.1CT 4 g. cm." 3 , where the radius rj up iter is 
the semi-major axis of the planets orbit. This density is above p\ and 
might be above pc- Taking a fluid to model the planets is an unusual step, 
but the alternative of seeking an n— body solution to the field equations is 
not viable because even the 2— body solution is not known. Looking at con- 
stant galactic rotation curves one might try an approximation. As noted in 
the last paragraph of section 5 of [ 14 1 : 

For constant circular velocities over a large distance it is neces- 
sary to have an approximately logarithmic potential. Thus the 
metric will have an approximately logarithmic term. The Rie- 
mann tensor is constructed from the second derivatives of the 
metric and the square of the first derivatives of the metric. For 
a logarithmic potential these will both be of the order r~ 2 and 
thus a linear analysis might not be appropriate. 

This suggests that only an approach using an exact solution will work. 
One can assume that the system under consideration can be modeled by 



a static spherically symmetric spacetime with line element |23[ Construct- 
ing the geodesies using Chandrasekhar's (1983) |37j method, the geodesic 
Lagrangian is given by 



2£ = -Ct 2 + Ar 2 + Bd 2 + B sin 2 d<j> z . (57) 
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The momenta are given by 

dC 

Pa = jr- (58) 
dx a 

and are 

p t = -Ci, p r = Ar, p e = BO, p^ = B sin 2 9<j). (59) 
Euler's equations are 

Pa = d a C (60) 

For static spacetimes with dtA = dtB = dtC = 0, giving Qjjr = so that the 

— n i- 

dr 



time component of the Euler equation 6C gives = 0, integrating 



Pt = C— = E a constant along each geodesic. (61) 
ar 



Similarly by spherical symmetry one can take d^A = d^B = d^C = 0, giving 

dp<p 
dr 



jJi = so that the <f> component of the Euler equation Bol gives = 0, 



integrating 

For the 9 component 



p<f, = Bsin 9— = a constant. (62) 



^ = B S in0cos0^, (63) 
o9 dr 



the Euler equation 6C is 



± pe = ^B6 = ^- = Bsin (9 cos A (64) 
dr dr 39 dr 

choosing to assign the value tt/2 to 9 when 9 is zero, then 9 will also be zero; 
and 9 will remain constant at the assigned value. The geodesic is described 
in an invariant plane which can be taken to be 9 = tt/2. Equation 62 now 
gives 

P0 = Bcj) = ia constant along each geodesic (65) 

where L is the angular momentum about an axis normal to the invariant 
plane. Substituting into the Lagrangian 

E 2 I? 
- — + Ar 2 + — = 2£ = -lor0, (66) 
U B 

where 2C = —1 or depending on whether time- like or null geodesies are 
being considered. Rearranging 

Af 2 = -— + — + 2£. (67) 
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Taking r to be a function of eft instead of r and using 



65 



gives 



dr\ 2 B B 2 E 2 \ , . 

now letting 

1 , s 

u = - 69 

r 



as in the usual Newtonian analysis 

seeking a thermodynamic interpretation one can substitute the enthalpy h 
for the lapse C = h 2 , but A and B are still arbitrary so that this is not 
pursued. Inserting the Kottler (Schwarzschild solution with cosmological 
constant) values of the metric 

B = r\ = I-l-*2 + ^, (71) 
A r 6 



and taking 2C = — 1 for time-like geodesies equation |70| becomes 

du\ 2 9 o 2mii l-£^ 2 A A . , 

— = -u 2 + 2mtf ri H ^ ^ (72) 

.64 J L 2 L 2 3u 2 L 2 3 V ' 

which is equation (4) of Reference [29|, the last term suggesting the pos- 
sibility of constant rotation curves. One can if investigate if there is any 
adjustment of the Newtonian potential which will give constant geodesies, 
as required for galactic rotation. Taking (c.f. Will (1993)|4l|l eq.4.6) 

ga ~- 1 + 2U, (73) 

where U is the Newtonian gravitational potential. Now in assume addition- 
ally the particular form for a spherically symmetric spacetime 

B = r 2 , A=I« T J^«l + 2^ (74) 
inserting in [7G and expanding for small U everywhere 



du\ 2 1-E 2 „. . 2 \J 

n 2 (l + 2C/) + -y (75) 



- 



■■>) L 2 v ' L 2 
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In particular one might expect that constant rotation is given by the middle 
term so that 

— u 2 (l + 2U) = a a constant, (76) 
rearranging for U we find 



U 



a 



(77) 



this suggests that the correct addition to U to produce constant curves is a 
function in r 2 , this is given by the addition of a cosmological constant and 
such a spacetime is given by Kottler's solution 71. One might ask what is 



the next simplest space-time after on with a cosmological constant which 
has an r 2 increasing potential and perhaps this is the interior Schwarzschild 
solution. This can be thought of as modeling the halo of a galaxy with the 
interior Schwarzschild solution and calculating the geodesies to see if they 
give constant motion. Newtonian modeling has been done by Binney and 
Tremaine (1987) ^2j. For the interior Schwarzschild solution Adler, Bazin, 
and Schiffer, (1975) jllj equation number 14.47 one has 



A 



C 



inserting into 70 one gets 




for r < ro, R 



3c 2 



(78) 



du\ 



1 1 

I^ + I 2 



-1 + 



i 2 R 2 



+ 



2E 



u 2 B 2 - 



•i 2 r 2 



u 2 R 2 -l 



V 



(79) 



The -Jj = ^rterm can be thought of as giving constant rotation curves 



proportional to the halo density v c oc p. What one expects from the Tully- 
Fisher (1977) j£| relationship is that oc L oc M. There might be an exact 
solution in Delgaty and Lake (1998) |^] which will model this more closely. 



19 



5.3 Rates of Decay. 

In general one can ask what exact solutions to gravitational field equations 
give what rate of decay. This problem could also be studied numerically. 
The rate of decay of scalar fields has been discussed in the last paragraph 



of the introduction of 12]. Including other fields one roughly gets the rates 
of decay: type B conformal scalars > perfect fluids > the gravitational field 
> type O scalars > electromagnetic fields > type A conformal scalars > 
coupled and interacting fields and fluids. The type B conformal scalars and 
perfect fluids are not usually asymptotically flat. Of course, for example, 
one would expect there to be conformal scalar solutions which are neither 
type A or B and these have unpredictable rates of decay, so this ordering is 
not absolute. 

5.4 The Spacetime of Elementary Particles. 

Rates of decay are not only important on long distance scales. As pointed 
out in the second paragraph of the introduction of the exact solution 
of the spherically symmetric spacetime of the Klein-Gordon-Einstein equa- 
tions is not known, except in the massless case where the static spherically 
symmetric field equations R a b = 14> a 4>b have the solution 

ds 2 = exp y ^ dt 2 — ^ exp ( J cosech 4 ^— ^ dr 2 

„2„„„ /2m\ 2 fv\ .2, 



-if exp [-1 cosech^ i d 2 S, = a/r (80) 

where rj 2 = m 2 + a 2 and m is interpreted as the mass and a the scalar 
charge, see for example, Roberts (1985) ||. Only the massless exact so- 
lution is known so that the exact modification of the shape of the Yukawa 
potential for the meson is not known: it proves difficult to approximate. The 
spacetime of mesons has been discussed by Fisher (1948) |p5| , Ross (1972) 
[4§, Nagy (1979) |47|], and Ho (1995) @. The Yukawa potential was in- 



vented, Yukawa (1935) ||49j, Landau (1990) |50fl , to account for the 7r-meson 
as the exchange quantum in the force between two nucleons; this is by anal- 
ogy with electromagnetism where it is the exchange of a photon that is the 
origin of the electric and magnetic forces between electrons. The exact form 
of the potential is V = — (1/r) exp(— r/m^) times a function which involves 
the relative spin orientations. The Yukawa potential is only an approxima- 
tion as Quantum Chromodynamics is really the theory of strong interaction; 
the 7r-meson or pion successfully describes the residual force, and is thought 
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to work up to momentum transfer of about 1 GeV/c. In strong interaction 
there is also evidence of a linear confining potential. The Yukawa potential 
for mesons can be measured using form factors, several mesons are needed, 
see for example Gross, Van Orden and Holinde (1992) |5l]]. The hypothetical 
Higgs particle also has a Yukawa potential. In order to measure the Higg's 
Yukawa potential one needs to measure the coefficients of the H 3 and H 4 
terms in 

V = M H2/2 H 2 + (M HV2v ) H 3 + (M H2/2v2 ) H\ (81) 

and verify that the coefficients are as expected. At present there are no 
measurements of these terms. To measure them one would need to observe 
multi-Higgs production: this is further discussed in Djouadi el al (1999) 
[p2| . In the standard model the Higgs coupling enters by a quartic coupling: 
however in supersymmetric theories the quartic couplings are connected to 
gauge couplings which are known, so that in supersymmetric models it is 
easier to calculate the coefficients. 

5.5 What Happens on Long Distance Scales. 

If asymptotic flatness is incorrect then what does happen on long scales? 
Non-asymptotic flatness introduces the problem of what happens as the 
potential goes to infinity - does it increase for ever? What could happen 
is that the one body problem becomes inappropriate: one needs a solution 
which takes into account two or more bodies. In particular for the solar 
system if there is a growing term in the potential one might take that it has 
stopped growing well before the next star. Torbett and Smoluchowski (1984) 
[f)3"f argue that there are bodies orbiting the Sun at 10 5 AU ~ 5 x 10 _3 pc, 
which might be a maximum orbiting distance. Puyoo and Jaffel (1998) 
study the interface between the heliopause and the interstellar medium, this 
is at about 10 3 AU and they find a high interstellar hydrogen density of 
0.24 ± 0.05 g.cm. -3 , a proton density of 0.043 db 0.005 g.cm. -3 , a helium 
density of (2.7 ± 0.5)10 -2 g.cm. -3 , and so forth. One consequence of these 
non-vanishing densities is that in gravitation, as in quantum field theory, it 
becomes difficult to say what a vacuum is and whether it has energy, see 
Roberts (2000) |55|. One can ask what sort of metric describes spacetime 
at various distances from the Sun, and it seems that some sort of onion 
model is called for. The standard picture is that the interior Schwarzschild 
solution is matched to the exterior Schwarzschild solution as in Adler et 
al (1975) §14.2., and then match the exterior Schwarzschild solution to a 
Friedman model with a specific equation of state as in Stephani (1985) |)(| 
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§27.3. Perhaps there should be more than three regions. The sun has a mean 
density of psun = 1.409g.cm. -3 Allen (1962) [p^| ; however its density varies 
considerably depending on its distance from the centre from Allen (1962) 
f37f table on page 163; there is a big jump at about half its radius, which 
can be modeled by a dense core, so perhaps two interior solutions are needed 
to describe it. Dziembowski et al (1990) and Basu et al (2000) |59|, use 
inversion techniques to show that the sun has many layers with different 
speeds of sound and densities. The solar system splits up into three regions: 
the inner where the general relativistic corrections to Newtonian theory are 
needed, the middle where Newtonian theory works, and the outer where a 
term explaining the irregularity in Pluto's orbit is needed. Next one needs 
a metric to describe the effect of local stars, then of the galaxy, and then of 
groups of galaxies, the Robertson- Walker cosmological region comes next, 
and after this perhaps a chaotic region. One can ask if a particle, say at 1 
parsec from the Sun is not in a flat region what is it that causes the most 
deviation from flatness. For simplicity assume that a Newtonian potential 
will give correct ratios between the contributions, so that the quantity <j)/G = 
M/R is calculated in units of the Sun's mass over parsecs. A parsec from 
the Sun is about as isolated as a particle in the nearby galaxy could be 
expected to be. The deviation from flatness of the metric is approximately 
given by equation 73 with U = 4>. The quantities in Allen (1962) p7[| 
§132,133,135,136, for the masses and distances associated with the local star 
system (Gould belt), the galaxy, the local group of galaxies, the Universe are 
used. Working to the nearest order of magnitude, the local star system has 
diameter 1,000 pc. and mass 1 x 10 8 Msum assuming the Sun is near the 
edge gives the potential M/R 10 5 Mg un pc. . The galaxy has diameter 
25 kpc. and mass 1.1 x 10 11 Mg un , but the distance of the Sun from the 
centre is 8.2±0.8 kpc, using this distance M/R m 10 8 M Sm pc. _1 . The local 
group of galaxies consists of 16 galaxies, suggesting an approximate mass of 
10 12 Msun, whose centre is 0.4 Mpc. away giving M/R « 10 7 Ms U nV c - 1 ■ 
Van den Berg (1999) |6(| finds 35 local group members and mass M^q = 
(2.3 ± 0.6) x 10 12 M Su n; and that the zero surface velocity, which separates 
the local group from the field that is expanding with the Hubble flow, has 
radius Ro = 1.18 ±0.15 Mpc The Universe has a characteristic length sale 
R = c/H ~ 3,000 Mpc. and the mass of the observable Universe is 10 54 g., 
again one can form a ratio M/R, but it has no direct meaning because 
of homogeneity, one finds M/R ~ 10 11 . To compare with the potential 
on the surface of the Earth note that the Earth's mean radius REarth ~ 
6xl0 3 Km. = 2xl0~ 12 pc. and has mass M Earth f» 6xl0 27 g. = 3xl0~ 6 M Su n, 
giving M/R ~ 10 6 Mg^pc. -1 for the contribution from the Earth's mass, 
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M/R ~ 10 5 MsunPC. 1 for the contribution from the Sun's mass. Collecting 
these results together gives the ratios 

10 6 : 10 5 : 1 : 10 5 : 10 8 : 10 7 : 10 11 (82) 

This suggests that either the Newtonian approximation is not appropriate, 
that asymptotic flatness is not a physical notion, or both. For the onion 
model it suggests that the metric describing the effect of local stars, the 
galaxy, and the local group of galaxies might not be needed because of 
the Universe's higher ratio. Another approach to what sort of notions are 
useful in describing stellar systems is as follows. A priori one would not 
wish to exclude the possibility that near the centre of the galaxy there are 
stellar systems consisting of several stars, many planets, many asteroids and 
comets, lots of dust, and which are close say only a light year away from 
other stellar systems. Dynamics for such a stellar system perhaps could 
still be calculated in some regions, but there are no notions of a one-body 
system, vacuum field equations, or asymptotic flatness to use in an explicit 
manner. 

5.6 Mach's Principle. 

Mach's principle can be formulated in many ways: Barbour and Pfister 
(1995) || p.530 list 21, Bondi and Samuel (1997) || list 10. Different 
formulations can lead to contradictory conclusions: for example, Bondi and 
Samuel's (1997) |62|] Mach3 and MachlO give rise to diametrically opposite 
predictions when applied to the Lense-Thiring effect. A Newtonian formu- 
lation has equations which can be used to describe dynamics rather than 
recourse to dark matter, Roberts (1985) Q. Lack of asymptotic flatness 
suggests that a system cannot be isolated. This is unlike thermodynamics 
where isolated heat baths are ubiquitous, and unlike electrostatics where 
the charge inside a charged cavity can be zero. So why should a Minkowski 
cavity in a Robertson- Walker universe be excluded? Field equations and 
junction conditions allow this to be done, it has to be excluded by principle. 
The answer is that it is different from electrostatics as gravitation is monopo- 
lar in nature. Any departure from homogeneity in the exterior region to a 
charged cavity would mean a change in charge which would quickly attract 
the opposite charge and cancel out: however in the gravitational case this 
does not happen, a change in homogeneity exterior to a Minkowski cavity (I 
think) would quickly change the spacetime from being flat. The above sug- 
gests a new formulation of Mach's principle: there are no flat regions 
of physical spacetime. What happens for an initial value formulation 
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of this is unclear: presumably it means that a well-defined initial surface 
does not develop into a surface part of which is flat. The above statement 
of Mach's principle is a particular case of the statement of Einstein (1953) 
|| , Ehlers (1995) ||, and Bondi and Samuel (1997) || Mach9 there 
are NO absolute elements: a fiat metric is an absolute element. 

5.7 Isolated Systems. 

Another way of looking at asymptotic flatness is to note that it implies that 
the solar system is isolated. Isolated systems seem to be an ideal which is 
appealed to in order to make problems soluble. The necessity of addressing 
soluble problems is discussed in Medewar (1982) H6q| . In practice an isolated 
system is only an approximation, there is always some interaction with the 
external world and for the assumption of an isolated system to work this 
must be negligible. The assumption that systems can be isolated appears 
through out science, but there appears to be no discussion of what this 
involves in texts in the philosophy of science. Three examples of isolated 
systems are now given. The first is photosynthesis: one can think of each 
leaf on a tree as an isolated entity with various chemical reactions happening 
independent of the external world, but this is only an approximation as the 
leaf exchanges chemicals with the rest of the tree so perhaps the tree should 
be thought of as the isolated system, further one can think of the entire 
biosphere as an isolated entity which converts 3 x 10 21 Joules per year into 
biomass from a total of 3 x 10 24 Joules per year of solar energy falling on 
the Earth, see for example Borisov (1979) §1.2.1 j67j. The second is in 
thermodynamics and statistical mechanics, here the isolatabilty of systems 
is taken as a primitive undefined concept, see for example Rosser (1982) ||5ql 
page 38. The third is of experiments where a single electron is taken to be 
isolated, Ekstrom and Wineland (1980) |j{J: the single electron is confined 
for weeks at a time in a "trap" formed out of electric and magnetic fields. 



6 The Tolman-Ehernfest Relation. 



6.1 The Radiation Fluid. 



For a radiation fluid 7 = §, and by [[?] the fluid index is 

UJ = (3p)4. 



(83) 



The Stefan-Boltzman law is 



P 




(84) 
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where T is the temperature and a is the Stefan-Boltzmann constant. Thus 



(85) 



Assuming the spacetime is static and admits the rotation free vector [5], 
equations and 85 give 



(86) 



thus showing that the lapse N is inversely proportional to the temperature 
T. This is the Tolman-Ehrenfest (1930 ) p3| ] relation. Lapse only spacetimes 
have been studied by Roberts (1994) [ffOjand Schmidt (1996) Q. For the 

gives 



non-static case 46 and 



Ng 



(3)1 



(87) 



7 The Geometric-thermodynamic equation and Cos- 
mic Censorship. 

7.1 Scalar Field Solutuions. 

It is known that spherically symmetric asymptotically flat solutions to the 
Einstein massless scalar field equations do not posses event horizons, both in 
the static case Roberts (1985) |J and in the non-static case Roberts (1996) 
|jl2|| . Massless scalar field solutions are equivalent to perfect fluid solutions 
with 7 = 2 and U a = (/> a (— C C ) _ 2 ; for the above scalar field solutions the 
vector field is not necessarily timelike so that the perfect fluid correspon- 
dence does not follow through. It can be argued that an asymptotically 
flat fluid would be a more realistic model of a collapsed object, because a 
fluid provides a better representation of the stress outside the object. In 
the spherically symmetric case a global coordinate system of the form ^ 
can be chosen and a necessary condition for there to be an event horizon is 
that, at a finite non-zero value of r, C — > oo. From |28|, |29|, 



and 38 it is 



apparent that this only occurs from some exceptional equations of state and 
values for the fluid density. Relaxing the requirement of spherical symmetry 
equations 55 and 56 show that for there to be a null surface iV — > 0, or 
u — ► oo; however the derivation of both ^ and 56 requires the vector 
and components of this diverge as N — > 0, also to show that |55| and |56| hold 



globally it is necessary to show that the coordinate system 39 can be set 
up globally. The above suggests that it is unlikely that spacetimes with a 
perfect fluid present have event horizons except in contrived circumstances. 
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